Synergistic Effect of P Doping and Mo-Ni-Based Heterostructure Electrocatalyst for Overall Water Splitting

Heterostructure construction and heteroatom doping are powerful strategies for enhancing the electrolytic efficiency of electrocatalysts for overall water splitting. Herein, we present a P-doped MoS2/Ni3S2 electrocatalyst on nickel foam (NF) prepared using a one-step hydrothermal method. The optimized P[0.9mM]-MoS2/Ni3S2@NF exhibits a cluster nanoflower-like morphology, which promotes the synergistic electrocatalytic effect of the heterostructures with abundant active centers, resulting in high catalytic activity for the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER) in alkaline electrolyte. The electrode exhibits low overpotentials and Tafel slopes for the HER and OER. In addition, the catalyst electrode used in a two-electrode system for overall water splitting requires an ultralow voltage of 1.42 V at 10 mA·cm−2 and shows no obvious increase in current within 35 h, indicating excellent stability. Therefore, the combination of P doping and the heterostructure suggests a novel path to formulate high-performance electrocatalysts for overall water splitting.


Introduction
The excessive use of traditional fossil fuels has caused a series of problems, such as global warming and environmental pollution, and the dependence on nonrenewable energy sources has led to resource depletion [1,2]. Therefore, the exploitation of fresh and clean energy sources is an efficacious strategy to accomplish low-carbon emissions and mitigate environmental problems [3]. Hydrogen is an ideal clean fuel for solving energy and environmental issues [4,5]. Currently, most sustainable hydrogen processes are not clean enough; nevertheless, the electrochemical water-splitting method can be used more prevalently as a new approach to produce high-purity H 2 . However, the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER) involve complex multielectron transfer mechanisms, resulting in a thermodynamic potential of 1.23 V. The high potential is considered a major challenge for water splitting [6]. Noble metal Pt-based [7] and Ru/Ir-based oxides have been considered ideal HER and OER electrocatalysts because of their good catalytic intrinsic activity [8]. However, their high price, scarcity, and inferior stability limit their commercial applications [9,10]. Thus, the production of inexpensive water-splitting electrocatalysts with high catalytic activity, conductivity, and durability is highly desirable.
Based on the construction of heterostructures, heteroatoms [27][28][29][30][31] are usually introduced to increase the number of active sites, optimize the electronic structure, and improve the reaction kinetics for the HER and OER. In recent years, V-doped Ni 3 S 2 nanowire arrays (V-Ni 3 S 2 ) have been prepared using a hydrothermal method [29], and the results showed that a large number of free carriers produced by V doping improved the catalytic performance of Ni 3 S 2 . Ma et al. [32] prepared N-doped hierarchical MoS x /Ni 3 S 2 -heterostructured nanowires on nickel foam (N-MoS x /Ni 3 S 2 -4@NF) for use as hydrogen evolution catalysts. They confirmed that N doping enhanced the conductivity and charge transfer rate of MoS x . P doping can expose more active sites and enhance the intrinsic conductivity and electrocatalytic activity of the catalysts [33,34]. Xue et al. [33] introduced P into the MoS 2 matrix using the solid-phase ion exchange method, and the synthesized P-MoS 2 catalyst exhibited good electrocatalytic performance for both the HER and OER. The P-MoS 2 /Mo 2 C electrocatalyst, which was prepared via high-temperature phosphating in a H 2 atmosphere, exhibited highly efficient HER activity [34]. P doping also increases the number of active sites and reduces the adsorption free energy. Wang et al. [35] used P-induced electron density modulation to activate the basal plane of WS 2 (P-WS 2 ), and the HER performance of the P-WS 2 nanowire arrays is superior to that of undoped WS 2 . Although the activity of these electrocatalysts has been vastly improved, P-doped electrocatalysts still require a lower overpotential for overall water splitting. Therefore, an effective approach to obtain highly efficient electrocatalysts is to further improve the synergistic effect of P doping and interface heterostructures, which could reduce the overpotential and increase the number of active sites for both the HER and OER processes.
Herein, a self-supporting P-doped MoS 2 /Ni 3 S 2 @NF heterostructure electrocatalyst was synthesized using a one-step hydrothermal method. Through the optimization and regulation of P doping, the synergistic effect of P doping and the interface heterostructure of the electrocatalyst resulted in remarkable electrocatalytic activity toward the overall water splitting. This study provides a new strategy for Mo-Ni-based electrocatalysts for efficient water electrolysis.

Experimental
2.1. Preparation of P [xmM] -MoS 2 /Ni 3 S 2 @NF and MoS 2 /Ni 3 S 2 @NF (NH 4 ) 2 MoS 4 (41 mg) and (NH 4 ) 2 HPO 4 (0.3, 0.6, 0.9, 1.2, and 1.5 mM) were used as Mo, S, and P sources, respectively, which were fully dissolved in 35 mL DMF solvent by magnetic stirring and ultrasound to form a homogeneous solution. The solution and the pretreated carrier NF (1 × 4 cm 2 ) were then transferred to a 50 mL polytetrafluoroethylene (Teflon)-lined stainless steel autoclave and heated at 200 • C for 12 h. Finally, the P [xmM] -MoS 2 /Ni 3 S 2 @NF catalyst electrodes with different initial P concentrations were obtained by repeatedly washing with deionized water and ethanol, and then drying at 60 • C for 10 h. A similar preparation method was used for MoS 2 /Ni 3 S 2 @NF, except that (NH 4 ) 2 HPO 4 was no longer added.

Preparation of MoS 2 @NF
When synthesizing MoS 2 powder, the mixed solution of (NH 4 ) 2 MoS 4 (41 mg) and DMF (35 mL) was directly transferred to a Teflon-lined autoclave. After heating at 200 • C for 12 h, the black powder at the bottom of the autoclave was collected, washed several times, centrifuged, and dried. The MoS 2 powder (7.13 mg) with the same loading as the target electrode was fully dissolved with ethanol, deionized water (V H2O /V ethanol = 5/4), and Nafion binder (10 µL) to form a homogeneous solution, which was uniformly coated on NF (1 × 1 cm 2 ) to prepare the MoS 2 @NF electrode.

Preparation of Ni 3 S 2 @NF
The preparation process is similar to that of P-MoS 2 /Ni 3 S 2 @NF. The NF was immersed in the homogeneous solution formed by CH 3 CSNH 2 (91 mg) and DMF (35 mL), and then heated in a Teflon-lined stainless steel autoclave at 200 • C for 12 h. Finally, the Ni 3 S 2 @NF electrode was prepared after cleaning and drying.

Phase and Structural Analyses of P-MoS 2 /Ni 3 S 2 @NF
A simple one-step hydrothermal process was used to prepare the P-MoS 2 /Ni 3 S 2 @NF catalyst for the overall water splitting by controlling the initial concentration of the P doping. Figure 1 shows a schematic of the synthesis process. The three-dimensional porous NF structure ( Figure S1g)  times, centrifuged, and dried. The MoS2 powder (7.13 mg) with the same loading as the target electrode was fully dissolved with ethanol, deionized water (VH2O/Vethanol = 5/4), and Nafion binder (10 µL) to form a homogeneous solution, which was uniformly coated on NF (1 × 1 cm 2 ) to prepare the MoS2@NF electrode.

Preparation of Ni3S2@NF
The preparation process is similar to that of P-MoS2/Ni3S2@NF. The NF was immersed in the homogeneous solution formed by CH3CSNH2 (91 mg) and DMF (35 mL), and then heated in a Teflon-lined stainless steel autoclave at 200 °C for 12 h. Finally, the Ni3S2@NF electrode was prepared after cleaning and drying.

Phase and Structural Analyses of P-MoS2/Ni3S2@NF
A simple one-step hydrothermal process was used to prepare the P-MoS2/Ni3S2@NF catalyst for the overall water splitting by controlling the initial concentration of the P doping. Figure 1 shows a schematic of the synthesis process. The three-dimensional porous NF structure ( Figure S1g) displayed a role as a catalyst carrier and provided the Ni source, owing to its good electrical conductivity. The P-doped MoS2/Ni3S2@NF nanoflower heterostructure was formed via a hydrothermal reaction with MoS4 2-and PO4 3-in (NH4)2MoS4 as well as (NH4)2HPO4, respectively. The surface morphology and structure of the catalysts were analyzed using scanning electron microscopy (SEM). The SEM image of the MoS2/Ni3S2@NF catalyst exhibits an irregular particle stacking structure (Figures 2a,b and S1f). When the heteroatom P was introduced, the morphology and structure changed, and the low-and high-resolution SEM images of the P[xmM]-MoS2/Ni3S2@NF (x = 0.3, 0.6, 0.9, 1.2, and 1.5 mM) doped with different P concentrations are shown in Figures 2c,d and S1a-e. When the initial P concentration was 0.3 mM, the granularity significantly decreased and lamellar growth was observed ( Figure S1a). A small amount of nanoflowers (approximately 50 nm) began to form at a P concentration of 0.6 mM ( Figure S1b). Finally, when the P concentration was increased to 0.9 mM, the nanoflowers were uniformly clustered to form a nanosphere heterostructure with a size of approximately 200 nm (Figures 2c,d and S1c) and the spherical surface was distinctively covered with lamellae. As the P concentration increased, the nanoflowers gradually disappeared ( Figure S1d) and finally became undulating cells ( Figure S1e) owing to the excessive attachment. The results demonstrate that the content of P doping can affect the specific surface area of the catalytic materials The surface morphology and structure of the catalysts were analyzed using scanning electron microscopy (SEM). The SEM image of the MoS 2 /Ni 3 S 2 @NF catalyst exhibits an irregular particle stacking structure (Figure 2a,b and Figure S1f). When the heteroatom P was introduced, the morphology and structure changed, and the low-and high-resolution SEM images of the P [xmM] -MoS 2 /Ni 3 S 2 @NF (x = 0.3, 0.6, 0.9, 1.2, and 1.5 mM) doped with different P concentrations are shown in Figure 2c,d and Figure S1a-e. When the initial P concentration was 0.3 mM, the granularity significantly decreased and lamellar growth was observed ( Figure S1a). A small amount of nanoflowers (approximately 50 nm) began to form at a P concentration of 0.6 mM ( Figure S1b). Finally, when the P concentration was increased to 0.9 mM, the nanoflowers were uniformly clustered to form a nanosphere heterostructure with a size of approximately 200 nm (Figure 2c,d and Figure S1c) and the spherical surface was distinctively covered with lamellae. As the P concentration increased, the nanoflowers gradually disappeared (Figure S1d) and finally became undulating cells (Figure S1e) owing to the excessive attachment. The results demonstrate that the content of P doping can affect the specific surface area of the catalytic materials [36], which can effectively modify the morphology and structure of Mo-Ni-based materials, while the high content of P doping changes the structure of the materials to a great extent [33]. In addition, the low-resolution SEM image in Figure S1g shows that the surface of the NF was initially smooth. However, the structure of the NF scaffold became rough after the P [xmM] -MoS 2 /Ni 3 S 2 synthesis (Figure S1a-f) due to an increase in the number of active attachment sites on the catalyst. The internal structural characteristics of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF were further analyzed using transmission electron microscopy (TEM). Figure 2e,f and Figure S2a,b show typical nanoflower structures and the presence of irregular flakes around the sphere, which may be due to the low crystallinity of the crystal [34]. The crystal structures of the nanoflowers were further explored using high-resolution TEM (HRTEM). The lattice fringes of the inner and outer lamellar structures of the spheres are shown in   [36], which can effectively modify the morphology and structure of Mo-Ni-based materials, while the high content of P doping changes the structure of the materials to a great extent [33]. In addition, the low-resolution SEM image in Figure S1g shows that the surface of the NF was initially smooth. However, the structure of the NF scaffold became rough after the P[xmM]-MoS2/Ni3S2 synthesis ( Figure    The crystal-phase structure of the prepared P [xmM] -MoS 2 /Ni 3 S 2 @NF catalyst was characterized using X-ray diffraction (XRD), and the results are shown in Figure 3a and Figure S4a. The XRD patterns of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF, MoS 2 /Ni 3 S 2 @NF, and Ni 3 S 2 @NF catalyst electrodes are shown in Figure 3a. The characteristic peaks of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF at 2θ = 21.  Figure S4b displays the XRD pattern of the MoS 2 powder. Owing to the low crystalline state [30,38], the MoS 2 diffraction peak is relatively broad, which is consistent with the HRTEM results. No obvious peak shift was observed between the Ni 3 S 2 and Ni after P doping, indicating that the introduction of P did not change the crystal structure of the MoS 2 /Ni 3 S 2 .

Electrocatalytic HER Activity
The HER performance of the P[0.9mM]-MoS2/Ni3S2@NF electrocatalysts was evaluated using a three-electrode system in 1 M KOH alkaline electrolyte, and the hydrogen evolution activities of MoS2@NF, Ni3S2@NF, MoS2/Ni3S2@NF, Pt/C@NF, and NF were tested for comparison. The obtained catalysts, a carbon rod, and Hg/HgO were used as the working, counter, and reference electrodes, respectively. Figure 4a displays the linear sweep voltammetry (LSV) curves of the various catalyst electrodes after iR correction. Other than the noble metal Pt/C@NF catalyst, which exhibited the best activity, the optimized P[0.9mM]-MoS2/Ni3S2@NF electrode also exhibited remarkable HER activity. When the same cathode current was reached, the applied driving voltage was the lowest and the current density increased sharply. Figure 4b exhibits the overpotential (η10) required for each electrode when the current density reached 10 mA·cm −2 . The P[0.9mM]-MoS2/Ni3S2@NF exhibited a low overpotential of 86 mV, which is smaller than that of MoS2@NF (139 mV), Ni3S2@NF (206 mV), MoS2/Ni3S2@NF (127 mV), and NF (225 mV). The substantial improvement in the HER activity indicates that the heteroatom P doping The chemical composition and valence of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst were further analyzed by X-ray photoelectron spectroscopy (XPS). The full spectrum displays the presence of elements such as Ni, Mo, S, and P (Figure 3b), which is consistent with the EDX results. The high-resolution XPS spectrum of Ni 2p is shown in Figure 3c. The P [0.9mM] -MoS 2 /Ni 3 S 2 @NF exhibits strong Ni 2p 3/2 and Ni 2p 1/2 peaks at binding energies of 856.31 and 874.18 eV, respectively. The energy difference of 17.87 eV between the two peaks indicates the coexistence of the Ni 2+ and Ni 3+ oxidation states [39], in which "Sat" is the oscillating satellite peak formed by the Ni-O bond [40]. Compared with the single-component Ni 3 S 2 @NF, the Ni 2p 3/2 and Ni 2p 1/2 peaks of the MoS 2 /Ni 3 S 2 @NF and P [0.9mM] -MoS 2 /Ni 3 S 2 @NF show a positive shift. Compared with MoS 2 /Ni 3 S 2 @NF, the prominent Ni 2p 3/2 and Ni 2p 1/2 peaks continue to shift positively by 0.23 and 0.28 eV after P doping. It was confirmed that the MoS 2 /Ni 3 S 2 heterostructure achieved interfacial charge transfer and distribution through the Mo-S-Ni bonds [41]. The high-resolution XPS spectrum of Mo 3d (Figure 3d) [42,43]. The signals at 229.08, 232.73, and 235.97 eV correspond to the characteristic peaks of Mo 3d 5/2 , Mo 3d 3/2 , and Mo 6+ of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF, respectively. Compared with the MoS 2 /Ni 3 S 2 @NF, the Mo 3d 3/2 and Mo 6+ peaks of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF show positive shifts of 0.41 and 0.52 eV, and the Mo 3d 5/2 peaks show negative shifts of 0.18 eV. Furthermore, the intensity of the peak of Mo 3d 5/2 is significantly enhanced, whereas that of Mo 6+ is weakened, which may be attributed to the reduction in the Mo valence due to the formation of P-Mo bonds by P doping [44]. In addition, a small S 2s peak was observed at 226.42 eV owing to oxidation, indicating the formation of S-O bonds. Figure 3e displays the high-resolution S 2p spectrum. The binding energies of the characteristic S 2p 3/2 and S 2p 1/2 peaks of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF were 162.23 and 162.97 eV, respectively, indicating the formation of the S 2− corresponding metal sulfides at the Ni-S and Mo-S bond sites [45]. Compared to the single-component Ni 3 S 2 @NF and MoS 2 , the corresponding peaks of the target electrode exhibit a certain energy shift. In particular, the characteristic S 2p 3/2 and S 2p 1/2 peaks of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF negatively shifted by 0.56 and 1.37 eV, respectively, compared to the MoS 2 /Ni 3 S 2 @NF. In addition, the intensity of the S 2p 1/2 peak is significantly enhanced, which corresponds to the results of the Mo 3d spectrum, indicating an increase in the Mo-S bond energy. In addition, a high-intensity peak of the residual SO 4 2− species (168.88 eV) was formed owing to S oxidation. Compared with MoS 2 /Ni 3 S 2 @NF, the positive shift of the Ni 2p, Mo 3d 3/2 , and Mo 6+ peaks and the negative shift of the Mo 3d 5/2 and S 2p peaks indicate that there is charge transfer between the MoS 2 and Ni 3 S 2 interface, which is beneficial to improve the electrocatalytic reaction process [42]. Figure 3f exhibits the high-resolution P 2p spectrum. The P [0.9mM] -MoS 2 /Ni 3 S 2 @NF spectrum exhibits strong peak signals at 133.33 and 134.27 eV, which correspond to P 2p 3/2 and P 2p 1/2 , respectively. This is possibly ascribed to the P-Mo and P-Ni bonds, which further verify the successful introduction of P into MoS 2 /Ni 3 S 2 @NF [33,46,47]. The XPS results show that there is a strong electronic interaction between the MoS 2 and Ni 3 S 2 species, suggesting that the establishment of the MoS 2 /Ni 3 S 2 heterointerface and P doping are beneficial for further enhancing the electronic interaction and synergistically optimizing the electronic structure. In addition, P can regulate the charge redistribution at the interface of the MoS 2 /Ni 3 S 2 heterostructure [34], thereby improving the charge transfer and reaction kinetics of Mo-S-Ni in electrocatalysis.

Electrocatalytic HER Activity
The HER performance of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF electrocatalysts was evaluated using a three-electrode system in 1 M KOH alkaline electrolyte, and the hydrogen evolution activities of MoS 2 @NF, Ni 3 S 2 @NF, MoS 2 /Ni 3 S 2 @NF, Pt/C@NF, and NF were tested for comparison. The obtained catalysts, a carbon rod, and Hg/HgO were used as the working, counter, and reference electrodes, respectively. Figure 4a displays the linear sweep voltammetry (LSV) curves of the various catalyst electrodes after iR correction. Other than the noble metal Pt/C@NF catalyst, which exhibited the best activity, the optimized P [0.9mM] -MoS 2 /Ni 3 S 2 @NF electrode also exhibited remarkable HER activity. When the same cathode current was reached, the applied driving voltage was the lowest and the current density increased sharply. Figure 4b exhibits the overpotential (η 10 ) required for each electrode when the current density reached 10 mA·cm −2 . The P [0.9mM] -MoS 2 /Ni 3 S 2 @NF exhibited a low overpotential of 86 mV, which is smaller than that of MoS 2 @NF (139 mV), Ni 3 S 2 @NF (206 mV), MoS 2 /Ni 3 S 2 @NF (127 mV), and NF (225 mV). The substantial improvement in the HER activity indicates that the heteroatom P doping promotes synergy between MoS 2 and Ni 3 S 2 . The Tafel slopes of the electrodes were calculated and fitted (Figure 4c). Compared with MoS 2 @NF (75.8 mV·dec −1 ), Ni 3 S 2 @NF (149.9 mV·dec −1 ), MoS 2 /Ni 3 S 2 @NF (47.7 mV·dec −1 ), and NF (214.8 mV·dec −1 ), the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF has the lowest Tafel slope of 28.5 mV·dec −1 , which is close to the Tafel slope of the Pt/C@NF (27.9 mV·dec −1 ), and the lowest onset potential (η onset ) of 61 mV. These results indicate that the kinetic activity of the HER is enhanced with an increase in the hydrogen production rate, which corresponds to the high-activity polarization curve shown in Figure 4a. The electrochemical impedance spectroscopy (EIS) curves in Figure 4d display the charge transfer resistance (R ct ) of the different catalysts. It can be seen from the semicircle diameter of the fitted Nyquist curve that the R ct of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF (0.73 Ω) is the smallest among the catalysts, except for Pt/C@NF. Thus, P doping is advantageous for accelerating electron transport on the NF surface for hydrogen evolution and improving the charge transfer kinetics in the HER, which correlates well with the electronic structure regulation measured by XPS [29,48]. The HER catalyst can be maintained with good stability, and Figure 4e displays the LSV curves of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF before and after 3000 cycles. When the driving current density reached 50 mA·cm −2 , the η 50 value of the polarization curve changed from 120 to 123 mV, the potential difference was small, and the curve was highly coincident. The electrode was also tested for long-term hydrogen evolution electrolysis at η = 115 mV for 50 h using chronoamperometry (Figure 4e). The stability curve of the electrolysis process exhibited a smooth trend, and the percentage change in the current density slightly decreased to 99.91%, indicating that the cathodic current loss was negligible. Based on the LSV results, the EIS curves of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF electrode before and after 3000 cycles were obtained. As shown in Figure 4f, the diameter of the fitted semicircle remained almost unchanged, with certain repeatability and coincidence. The smooth step obtained from the multi-current ( Figure S7a) and multipotential steps ( Figure S7b) tests indicates that the catalyst has good stability and durability for the electrocatalytic HER and can meet the practical application requirements of the industry. The effects of different P-doping concentrations of the P [xmM] -MoS 2 /Ni 3 S 2 @NF catalysts on the HER performance were also explored. Figure S5 displays a comparison of the HER performance of the P [0.3mM] -MoS 2 /Ni 3 S 2 @NF, P [0.6mM] -MoS 2 /Ni 3 S 2 @NF, P [0.9mM] -MoS 2 /Ni 3 S 2 @NF, P [1.2mM] -MoS 2 /Ni 3 S 2 @NF, and P [1.5mM] -MoS 2 /Ni 3 S 2 @NF. After optimization, the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst, with a P-doped amount of 0.9 mM, exhibited the best electrocatalytic HER activity. Based on these results, P doping can improve the activity of MoS 2 /Ni 3 S 2 @NF, and too little or too much P doping can slightly decrease the activity. mA·cm −2 , the η50 value of the polarization curve changed from 120 to 123 mV, the potential difference was small, and the curve was highly coincident. The electrode was also tested for long-term hydrogen evolution electrolysis at η = 115 mV for 50 h using chronoamperometry (Figure 4e). The stability curve of the electrolysis process exhibited a smooth trend, and the percentage change in the current density slightly decreased to 99.91%, indicating that the cathodic current loss was negligible. Based on the LSV results, the EIS curves of the P[0.9mM]-MoS2/Ni3S2@NF electrode before and after 3000 cycles were obtained. As shown in Figure 4f, the diameter of the fitted semicircle remained almost unchanged, with certain repeatability and coincidence. The smooth step obtained from the multi-current ( Figure S7a) and multi-potential steps ( Figure S7b) tests indicates that the catalyst has good stability and durability for the electrocatalytic HER and can meet the practical application requirements of the industry. The effects of different P-doping concentrations of the P[xmM]-MoS2/Ni3S2@NF catalysts on the HER performance were also explored. Figure S5 displays

Electrocatalytic OER Activity
The OER performance of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF electrocatalyst was evaluated under identical conditions. Firstly, the P-doping concentration was optimized, and the electrocatalytic OER performances of the P [0.3mM] -MoS 2 /Ni 3 S 2 @NF, P [0.6mM] -MoS 2 /Ni 3 S 2 @NF, P [0.9mM] -MoS 2 /Ni 3 S 2 @NF, P [1.2mM] -MoS 2 /Ni 3 S 2 @NF, and P [1.5mM] -MoS 2 /Ni 3 S 2 @NF were evaluated ( Figure S6). It was observed that the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF exhibited good OER activity, which is consistent with the HER activity. The OER performances of the MoS 2 /Ni 3 S 2 @NF, MoS 2 @NF, Ni 3 S 2 @NF, IrO 2 @NF, and NF were further investigated, and the results are shown in Figure 5. Figure 5a displays the LSV curves of the catalyst electrodes. The P [0.9mM] -MoS 2 /Ni 3 S 2 @NF had the lowest η onset , and its anode current increased sharply with increasing driving voltage, even exceeding the performance of the noble metal IrO 2 @NF catalyst. Figure 5b displays the overpotential of each catalyst at a current density of 100 mA·cm −2 (η 100 ). Compared with MoS 2 @NF (416 mV), Ni 3 S 2 @NF (439 mV), MoS 2 /Ni 3 S 2 @NF (347 mV), IrO 2 @NF (400 mV), and NF (597 mV), the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF had the lowest overpotential of η 100 = 312 mV. The Tafel slope of each catalyst is shown in Figure 5c. Compared with MoS 2 @NF (188.8 mV·dec −1 ), Ni 3 S 2 @NF (206.2 mV·dec −1 ), MoS 2 /Ni 3 S 2 @NF (129.2 mV·dec −1 ), and NF (261.4 mV·dec −1 ), the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF had the smallest Tafel slope of 47.2 mV·dec −1 . The EIS fitting curve in Figure 5d displays the R ct values of the different electrodes. The semicircle diameter of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF was the smallest, affording a minimum impedance resistance of R ct = 0.68 Ω. These results confirmed that P doping and MoS 2 /Ni 3 S 2 heterostructures synergistically improve the charge transfer rate in the electrolyte [49]. The OER stability of the optimized P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst was also evaluated. The LSV (Figure 5e) and EIS (Figure 5f) curves before and after 3000 cycles were coincident, the η 100 of the polarization curve increased from 312 to 318 mV at a current density of 100 mA·cm −2 , and the potential difference was small. The oxygen evolution electrolysis test performed at η = 290 mV for 50 h using chronoamperometry is illustrated in Figure 5e. The stability curve of the electrolysis process exhibits a smooth and gradual downward trend; the percentage change in the current density decreases to 99.76%, whereas the anode current loss is negligible. The multi-current ( Figure S7c) and multi-potential steps ( Figure S7d) over a short period of time were also examined. It can be seen that the steps formed by the step current and potential are smooth and stable, indicating that the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF has good stability and durability for the electrocatalytic OER. Table S1 summarizes and compares the HER and OER performances of the previously reported electrocatalysts. The performance of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF is better than that of most non-noble metal catalysts reported to date and it also exhibits good electrocatalytic activity for the HER and OER.
To further explore the effect of P doping on the electrocatalytic active sites of the MoS 2 /Ni 3 S 2 @NF, the electrochemical active surface area (ECSA) test was conducted. The scanning speed of 20-100 mV/s and current density of j 0.504V = (j a − j b )/2 were used as the horizontal and vertical coordinates, respectively. The double-layer capacitance (C dl ) values of the P [xmM] -MoS 2 /Ni 3 S 2 @NF, MoS 2 /Ni 3 S 2 @NF, MoS 2 @NF, Ni 3 S 2 @NF, and NF were obtained via calculations and fitting, and the results are shown in Figure S8 and Figure 6. Compared with the morphological and structural analyses by SEM and TEM, a large number of uniform flower-like spheres were formed on the surface of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst, which increased the number of exposed interfacial active sites and enlarged the active area. Therefore, the effective ECSA value obtained at the same scanning speed and potential range was larger, and the image was rounder and more symmetrical. The ECSA and C dl curves of the P [xmM] -MoS 2 /Ni 3 S 2 @NF ( Figure  S8) demonstrate that as the P-doping concentration increases, the C dl value of the P [xmM] -MoS 2 /Ni 3 S 2 @NF initially decreases and then increases to the maximum value of 122.5 mF·cm −2 (x = 0.9 mM), after which it decreases to 20.2 mF·cm −2 (x = 1.5 mM). The ECSA curves (Figure 6a-e) and C dl values (Figure 6f) of the different catalysts show that the optimized P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst has the largest active surface area. Compared with the C dl values of the MoS 2 /Ni 3 S 2 @NF (68.7 mF·cm −2 ), MoS 2 @NF (27.9 mF·cm −2 ), Ni 3 S 2 @NF (15.0 mF·cm −2 ), and NF (6.5 mF·cm −2 ), the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF exhibits a maximum C dl value of 122.5 mF·cm −2 . The increase in the C dl value indicates that the proper doping of heteroatom P can promote the uniform integration of MoS 2 and Ni 3 S 2 , enriching a large number of active sites between the interfaces of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst, which increases the effective active area for electrocatalytic hydrogen production [50,51] and improves the electrocatalytic activity. polarization curve increased from 312 to 318 mV at a current density of 100 mA·cm −2 , and the potential difference was small. The oxygen evolution electrolysis test performed at η = 290 mV for 50 h using chronoamperometry is illustrated in Figure 5e. The stability curve of the electrolysis process exhibits a smooth and gradual downward trend; the percentage change in the current density decreases to 99.76%, whereas the anode current loss is negligible. The multi-current ( Figure S7c) and multi-potential steps ( Figure S7d) over a short period of time were also examined. It can be seen that the steps formed by the step current and potential are smooth and stable, indicating that the P[0.9mM]-MoS2/Ni3S2@NF has good stability and durability for the electrocatalytic OER. Table S1 summarizes and compares the HER and OER performances of the previously reported electrocatalysts. The performance of the P[0.9mM]-MoS2/Ni3S2@NF is better than that of most non-noble metal catalysts reported to date and it also exhibits good electrocatalytic activity for the HER and OER.   In summary, the P[0.9mM]-MoS2/Ni3S2@NF exhibited good electrocatalytic performance for the following reasons. The electronegativity of a P atom is less than that of a S atom, so P doping can produce a positive charge in a sulfide structure, and the lone pair electrons In summary, the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF exhibited good electrocatalytic performance for the following reasons. The electronegativity of a P atom is less than that of a S atom, so P doping can produce a positive charge in a sulfide structure, and the lone pair electrons in a P atom can increase the electron spin density [52] and increase the active sites in the reaction process, thereby promoting the progress of the HER and OER. In addition, a TEM analysis shows that P doping can expose more edge structure defects, and the porous nanoflower structure forms a large specific surface area, which provides a sufficient reaction site for the HER and OER. It is worth noting that an SEM analysis shows that the introduction of a high concentration of P destroys the nanoflower structure of sulfides, thus reducing the electrical conductivity of Mo-Ni-based materials. Therefore, proper P content doping (2.06 wt%) can significantly improve the electrocatalytic performance of MoS 2 /Ni 3 S 2 . In addition, we also explored the electrocatalytic HER performance of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF in neutral (0.5 M PBS, pH = 7) and acidic (0.5 M H 2 SO 4 , pH = 0) electrolytes ( Figure S9). Unfortunately, the performance of its electrocatalytic OER is poor, so further research is still needed.

Overall Water-Splitting Performance
Based on the evaluation results of the electrocatalytic HER and OER ( Figure S10 and Table S2), the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst exhibits excellent activity for hydrogen and oxygen evolution. Thus, its performance as a bifunctional catalyst for overall water splitting in a two-electrode system of 1 M KOH alkaline electrolyte (Figure 7a) was examined. In the experiment, the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF||P [0.9mM] -MoS 2 /Ni 3 S 2 @NF coupling electrode was used as the anode and cathode of the electrolysis cell. During the watersplitting process with an applied driving voltage, large amounts of H 2 and O 2 bubbles were continuously generated at the cathode and anode of the device diagram, respectively (Figure 7c). As shown in Figure 7b, the overall water-splitting polarization curve of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF reached 10 mA·cm −2 at a driving voltage of 1.42 V, and the combined overpotential of the electrocatalytic water-splitting process was only 190 mV, which is significantly lower than the voltage required for the MoS 2 /Ni 3 S 2 @NF electrolysis cell (1.51 V). The considerable decrease in the driving voltage indicates that the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst promotes the reaction kinetics of the HER and OER [9]. The performance of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst is superior to that of most nonnoble metal bifunctional catalysts previously reported (Figure 7d), including MoS 2 -Ni 3 S 2 HNRs/NF (1.50 V) [25], Fe-MoS 2 /Ni 3 S 2 /NF-2 (1.61 V) [30], Ni 3 S 2 /MnS-O (1.54 V) [50], and as-anodic CoS x /Co (1.64 V) [53]. Figure 7c exhibits the chronoamperometry curve of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF operating for 35 h at 1.5 V. The results indicate that the electrolysis cell provides a stable current. During the long-term water-splitting process, the percentage change in the cathode and anode current densities was 98.9%, indicating that the water-splitting performance decreased slightly by 1.1%. The morphology of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF after overall water splitting was also examined using SEM ( Figure S11a,b), showing that most of the nanoflower structures were retained. The EDX elemental mapping spectrum (Figure S11c-g) confirmed the uniform distribution of Ni, Mo, S, and P on the catalyst surface. In conclusion, the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF exhibits excellent overall water-splitting performance and stability in alkaline solution.
In order to explore the mechanism of P doping in the electrocatalytic reaction process, the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF was characterized by XPS after a long-term OER test (Figure 8). The XPS full spectrum after the test correlates well with that before the test ( Figure  S12). Compared to the high-resolution spectrum of the initial P [0.9mM] -MoS 2 /Ni 3 S 2 @NF sample, the corresponding characteristic peaks of the Ni 2p (Figure 8a) and Mo 3d (Figure 8b) spectra shifted to the low-energy direction and became electron-rich centers. In contrast, the corresponding characteristic peaks of the S 2p (Figure 8c) and P 2p (Figure 8d) spectra shifted to the high-energy direction and became electron-poor centers. In Figure 8c, the intensity of the S 2p 3/2 and S 2p 1/2 peaks was significantly weakened, which may be attributed to the leaching of S during the electrocatalytic process. Compared with the spectrum before the test, the characteristic P 2p 3/2 (133.75 eV) and P 2p 1/2 (134.61 eV) peaks of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF after testing were positively shifted by 0.42 and 0.34 eV, respectively (Figure 8d). The above results show that the S atom with weaker electronegativity has a charge exchange with Ni and Mo with stronger electronegativity, which optimizes the internal electronic structure of MoS 2 and Ni 3 S 2 . At the same time, the P atom with the weakest electronegativity also has electron transfer, and the electrons transfer from the P atom to Ni, Mo, and S, so that P has more positive charge, resulting in an increase in the OER active center material [54]. It is confirmed that P doping is beneficial to promote the synergistic effect of MoS 2 and Ni 3 S 2 in the OER process. In summary, the optimized electronic configuration due to P doping contributed to the enhanced electrocatalytic performance of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF heterostructure.
the P[0.9mM]-MoS2/Ni3S2@NF catalyst promotes the reaction kinetics of the HER and OER [9]. The performance of the P[0.9mM]-MoS2/Ni3S2@NF catalyst is superior to that of most nonnoble metal bifunctional catalysts previously reported (Figure 7d), including MoS2-Ni3S2 HNRs/NF (1.50 V) [25], Fe-MoS2/Ni3S2/NF-2 (1.61 V) [30], Ni3S2/MnS-O (1.54 V) [50], and as-anodic CoSx/Co (1.64 V) [53]. Figure 7c exhibits the chronoamperometry curve of the P[0.9mM]-MoS2/Ni3S2@NF operating for 35 h at 1.5 V. The results indicate that the electrolysis cell provides a stable current. During the long-term water-splitting process, the percentage change in the cathode and anode current densities was 98.9%, indicating that the water-splitting performance decreased slightly by 1.1%. The morphology of the P[0.9mM]-MoS2/Ni3S2@NF after overall water splitting was also examined using SEM ( Figure  S11a,b), showing that most of the nanoflower structures were retained. The EDX elemental mapping spectrum (Figure S11c-g) confirmed the uniform distribution of Ni, Mo, S, and P on the catalyst surface. In conclusion, the P[0.9mM]-MoS2/Ni3S2@NF exhibits excellent overall water-splitting performance and stability in alkaline solution.  with the weakest electronegativity also has electron transfer, and the electrons transfer from the P atom to Ni, Mo, and S, so that P has more positive charge, resulting in an increase in the OER active center material [54]. It is confirmed that P doping is beneficial to promote the synergistic effect of MoS2 and Ni3S2 in the OER process. In summary, the optimized electronic configuration due to P doping contributed to the enhanced electrocatalytic performance of the P[0.9mM]-MoS2/Ni3S2@NF heterostructure.

Conclusions
In summary, a P-doped MoS 2 /Ni 3 S 2 @NF heterostructure electrocatalyst was synthesized using a facile one-step hydrothermal method. The results of various tests and characterizations showed that the nanoflowers morphology formed by the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF exposed more active sites, resulting in a large specific surface area and abundant active centers. In addition, the introduction of P promoted the synergistic effect of the MoS 2 /Ni 3 S 2 heterostructure interface, optimized the electronic structure, and effectively improved the transport speed of the carriers in alkaline media. The combination of these engineering system properties synergistically enhanced the electrocatalytic performance and reaction kinetics of the P [0.9mM] -MoS 2 /Ni 3 S 2 @NF catalyst for the HER, OER, and overall water splitting. When used as an electrode in an electrolysis cell, P [0.9mM] -MoS 2 /Ni 3 S 2 @NF only requires an ultralow driving voltage of 1.42 V for water splitting at 10 mA·cm −2 , which is superior to that of most bifunctional electrocatalysts reported to date. This study opens a novel avenue for designing highly active, stable, and economical multi-component P-doped heterostructure electrocatalysts.